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Abstract 

We use a relativistic mean field model to study the charge density distributions of 
exotic oxygen isotopes. Nonlinear isoscalar-isovector terms which are not constrained 
by the present data are considered and the pNN tensor coupling that does not affect 
the symmetry energy in the mean field model is included. Strong correlations between 
the neutron radius of 208 Pb and the charge radius of 12 > 13 > 23 > 24 are found. The pNN 
tensor coupling explicitly enhances the radius correlations between 23 > 24 and 208 Pb. 
This enhancement is due to the neutron occupation of the orbital 2si/a. The charge 
radius is sensitive to the change of the density dependence of the symmetry energy as the 
isotope goes towards the proton drip line. 
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The nuclear equation of state (EOS) is fundamentally important to describe the nuclear 
phenomenology. The EOS comprises of two components: the saturation property of sym- 
metric matter and the asymmetric energy. There exists a Cooster band for the saturation 
property of symmetric matter. However, the density dependence of the symmetric energy is 
much more poorly known to date[l, 2]. Therefore, the understanding of amount of properties 
of neutron stars, nuclei far from the (3 stability, and radioactive ion collisions exists theoretical 
uncertainties. Many factors, for instance, the isoscalar-isovector coupling terms[3, 4, 5, 6, 7], 
the isovector-scalar mesons and density dependent coupling constants [8, 9], can modify the 
behavior of the density dependence of the symmetry energy. Recently the density depen- 
dence of the symmetry energy has been extensively explored through the inclusion of the 
isoscalar-isovector coupling terms. 

The inclusion of the isoscalar-isovector coupling terms allows one modify the neutron 
skin of heavy nuclei without changing the saturation property of symmetric matter and a 
variety of ground-state properties that are well constrained experimentally [3]. Theoretical 
estimates may give an uncertainty of the neutron radius of 208 Pb by about 0.3fm. Recently, 
data to data correlations between the neutron radius of 208 Pb and neutron star properties 
have been established [3, 4, 5, 6]. Such correlations are also found in the binding energy of 
valence neutrons for some neutron rich nuclei [7]. Accordingly, an accurate measurement of 
the neutron radius for 208 Pb at the Jefferson Laboratory[10] that promises a 1% accuracy is 
very significant. 

With the development of the radioactive beams properties of many exotic nuclei have been 
exploredfll]. Many neutron and proton halo nuclei are identified through the measurement 
of the reaction cross sections and longitudinal momentum distributions of nucleus breakup. 
To obtain more definite information of halo nuclei, a pure probe that excludes the strong 
interaction is necessary. At RIKEN[12] and GSI[13], new facilities that provide electron 
scattering on exotic nuclei are under construction. These facilities will provide an elaborate 
explorations on the charge density distributions of proton halo nuclei. Eventually, the charge 
density distributions of all exotic nuclei are yet to be surely measured. Recently, Wang et.al. 
investigated elastic electron scattering on light exotic proton-rich nuclei theoretically [14]. 
On the other hand, the theoretical charge density distributions for these nuclei may exist 
uncertainties since the properties of light nuclei far from (3 stability are partly related to the 
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poorly known density dependence of the symmetry energy. 

The charge radius of heavy nuclei such as 208 Pb is little modified with the change of 
the density dependence of the symmetry energy although the neutron radius is apparently 
modified [3]. Due to the strong couplings of the neighboring neutron and proton shells in 
light nuclei, it is not so simple to say that the charge radius is still unchanged. In turn, as 
the proton skin emerges for the light proton-rich isotopes, we even conjecture that a relation 
between the neutron radius of 208 Pb and the charge radius of the proton-rich isotope (for 
instance, 12 0) will be established. It is aim of this study to constrain the density dependence 
of the symmetry energy through the investigation of the charge density distribution of the 
light nuclei far from the [3 stability. 

The relativistic mean-field (RMF) theory (for reviews, see[15, 16]) is powerful in describing 
the ground-state properties of nuclei from the proton drip line to the neutron drip line, since 
it can give dynamically rise to the spin-orbit interaction associated with shell effects. In this 
work, we take into account the isoscalar-isovector terms that modify the density dependence 
of the symmetry energy in the framework of the RMF models. In a RMF model, the pNN 
tensor coupling that does not change any properties of nuclear matter plays its role in finite 
nuclei. The interest to include the pNN tensor coupling is arisen by the following factors. On 
one hand, the isoscalar-isovector terms give modifications to the source term of the p meson 
(i.e., the related density), on the other hand, the tensor coupling is relevant to the gradient 
of the related density. We will investigate oxygen isotopes for their simple spherical shapes. 

The effective Lagrangian density is given as follows 

+\(d^a - m 2 a a 2 ) - -F^ + \ml^ 

where ip, a, ui, and bo are the fields of the nucleon, scalar, vector, and charge-neutral isovector- 
vector mesons, with their masses Mn ^m a ^m w ^ and m p , respectively. A p is the field of the 
photon. gi(i = a,co,p) and f p are the corresponding meson-nucleon couplings, r and T3 is 
the isospin Pauli matrix and its third component, respectively. F pu , B^, and are the 
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strength tensors of the u>, p mesons and the photon, respectively 

Fpv = d^u v - dvU^, = d^b 0l/ - d u b 0lJL , A^ u = d^A u - d v A^ (2) 
The self-interacting terms of a, w-mesons and the isoscalar-isovector ones are in the following 

-4g 2 p (A s g 2 a a 2 + A v g 2 u;^)b X (3) 

Since detailed procedures are easily found in references (e.g., see Ref.[16]) we only write out 
the radial Dirac equations and the tensor potentials from the pNN tensor coupling explicitly 

^ = {M * N (r) + V(r) + E a + Uar))F a -^-U$(r)-U§(r))G a (4) 
^ -- (M* N (r)-V(r)-E a -U?(r))G a + (--U%(r)-Ug(r))F a (5) 



dr 



and 



f f°° 

U * {T) = ~2M^ 9p r dnrl dGf> dr n) (P3(ri) " ^s 9 W + K 9 Wo)b Q ) (6) 

f2 



ttP( ^ f P f° a 2 dG p (r, n ) 



where 



PT(r) = £ 2 ^t±(G a (r)F a (r)) (7) 

a 

M* N (r) = M N + V a {r) 

V(r) = gUr)Mr) + 9 P bo(r)t + e(^)A (8) 
V a (r) = g a {r)a(r) 

with p3 the difference of the proton and neutron densities, t=l for proton and -1 for neutron, 
and G p (r,r\) the Green function of the p meson. The various density definitions follow the 
standard procedure [16]. The charge density is obtained by folding the proton density and 
the proton charge distribution[17] 

Pc (r) = J d 3 r'^- exp(- M |r - r'|) Pp (r'), p = (O.n^GeV (9) 

with p p the proton density. 



The pairing correlation is involved in non-magic nuclei using Bardeen-Cooper-Schrief 
(BCS) theory. We use constant pairing gaps which are obtained from the prescription of 
Moller and Nix [18]: A n = 4.8/iV 1 / 3 , A p = 4.8/Z 1 / 3 with N and Z the neutron and proton 
numbers, respectively. For odd nuclei, the Pauli blocking is considered. 

We perform calculations with the NL3 parameter set [19] and distinguish two cases with 
and without the pNN tensor coupling for various isoscalar-isovector coupling A„'s. For sim- 
plicity, the isoscalar-isovector coupling A s is set as zero. The symmetry energy at saturation 
density is not well constrained, and hence some average of the symmetry energy at satura- 
tion density and the surface energy is constrained by the binding energy of nuclei. For a 
given coupling A„, we follow Ref.[3, 7] to readjust the pNN coupling constant g p so as to 
keep an average symmetry energy fixed as 25.68 at /cp = 1.15 fm -1 . In doing so, it was 
found in Ref.[3] that the binding energy of 208 Pb is nearly unchanged for various A„'s. In 
our calculation, the total binding energy of light exotic nuclei with different A^'s is modified 
by the same order of 208 Pb. The pNN tensor coupling f p is taken as f p /g p = 6.1 which is 
used in the Bonn potentials [20] and is within 8.0 ±2.0 obtained from the light cone sum rules 
by Zhu[21]. Note that this ratio of f p to g p , coupled with the relatively large values of g p 
in many relativistic mean field models that are usually fixed by the empirical values of the 
symmetry energy, could overestimate the contribution of the rho tensor coupling, and that 
will be analyzed below. 

Table 1: Results for 24 with the NL3 parameter set for f p = and f p = 6.1g p , respectively. 
The total binding energy (Eb), the proton radius (r p ), and the neutron thickness (r p — r n ) 
are given for various A„'s. The experimental binding energy is 168.48MeV. 
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0.010 


169.74 


2.637 
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0.015 
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0.020 


168.93 
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0.626 


170.82 


2.655 


0.580 


0.025 


168.32 


2.646 


0.621 


170.21 


2.660 


0.569 


0.030 


167.72 


2.649 


0.616 


169.43 


2.667 


0.557 



At the vicinity of the line N = Z, the p meson together with the isoscalar-isovector 
term has almost negligible contribution. Consequently, we study the isoscalar-isovector term 
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contribution in nuclei close to the drip lines of the oxygen. Table 1 gives the binding energy, 
proton radius, neutron thickness for 24 with and without the pNN tensor coupling as A„ 
is changed. The binding energy of 24 shows a small reduction as the A„ is increased, while 
for 208 Pb the binding energy shows a weak increase [3]. The neutron thickness is significantly 
reduced with the increasing A„. Without the tensor coupling, the proton radius shows a 
much less modification compared to the neutron radius, which is similar to that for 208 Pb[7]. 
As the pNN tensor coupling is included, an apparent difference appears: the proton radius 
modification due to the A^ is largely enhanced, and almost half the contribution to the 
decrease of the neutron thickness of 23,24 O comes from the increase of the proton radius. We 
notice that the inclusion of the tensor coupling in 208 Pb just brings out negligible changes 
due to the quite small tensor potentials. 
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Figure 1: The correlation of the neutron thickness of Pb and that of ' O. 

A data to data correlation between the neutron thickness of 208 Pb and that of 23 > 24 
is displayed in Fig.l. The correlation is linear, and the pNN tensor coupling increases the 
correlation slope strongly. The relativistic models with a softer symmetry energy have a larger 
symmetry energy at low densities. It was pointed out in Ref.[7] that the softer models are the 
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first ones to drip neutrons that show a diffusive distribution with the larger valence neutron 
radius. A similar result for 23 ' 24 is obtained that the radius of the valence 2s 1 / 2 neutron 
increases with softening the symmetry energy. On the other hand, the change of the neutron 
radius is mainly from the modification of core neutron distribution instead of the diffusive 
neutron distribution that is often called as the neutron skins or halos. Consequently, owing to 
the coupling between the core neutrons and protons at the neighboring shells in light isotopes 
the modification of the neutron density distribution leads to a similar modification for the 
proton (charge) density distribution. In Fig. 2, the charge density for 24 is displayed. As 
seen in Fig. 2, the modification of the charge density distribution occurs explicitly in the core 
region with various A„'s, and the modification is enhanced due to the inclusion of the pNN 
tensor coupling. Further, we make a brief analysis of the overestimation of the contribution 
of the pNN tensor coupling due to the relatively large values of g p in many relativistic mean 
field models. As we decrease the value of f p by 15%, the slope of the neutron thickness 
correlation, shown in Fig.l, is just reduced by 8%. In this sense, the importance of the pNN 
tensor coupling in present results is not much overestimated. 
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Figure 2: The charge density distributions for 24 O with various cases. 

The explicit enhancement of the correlation for 23 > 24 as shown in Fig.l due to the pNN 
tensor coupling does not occur for isotopes with more neutrons. At the same time, with 
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Figure 3: The neutron spin-orbit potential for 22 > 23 . 24 0. The neutron occupation probability 
in the orbital 2si/ 2 , just slightly influenced by the A„ and f p , is 0.35, 0.5, 0.83 for 22 > 23 > 24 0, 
respectively. The lower one of the same type curves in each panel represents for the case 
fp/gp = 6.1, while the upper one is for f p /g p = 0. 



the reduction of the neutron number the contribution from the pNN tensor coupling drops. 
These phenomena can be explained by the contribution from the orbital 2s 1 / 2 - The orbital 
2si/2 generates large derivatives of the density distribution with respect to the radius near 
the node at the core region. A large tensor coupling contribution, corresponding to large 
tensor potentials, may come from the large derivatives appearing in integrals in Eq.(6). For 
isotopes with N > 16, the 2s x / 2 neutron distribution is squeezed due to the deeper binding, 
and a prominent cancellation occurs in the integrals over the density derivatives from the 
regions inside and outside the node, leading to dropping tensor potentials. The derivatives 
outside the node are small for 24 due to the extensive distribution, and the cancellation is 
slight so that large tensor potentials are produced. For isotopes with N < 16, the drop of 
the contribution from the pNN tensor coupling is related to a smaller occupation probability 
in 1s\j2- Fig. 3 displays the neutron spin-orbit potential XJ\ S that is obtained in an equation 
for the big component G a of the Dirac spinor by eliminating the small component F a from 
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Eqs.(4,5). The spin-orbit potential Ui s varies sensitively with the different 2s 1 / 2 occupation. 
If the Pauli blocking is not considered in 23 0, the 2s 1 / 2 occupation probability shows a small 
increment by about 0.1 that may produce an obvious move towards the spin-orbit potential 
for 24 0. As the occupation probability in 2si/ 2 rises, the difference arisen by the pNN tensor 
coupling goes up, and that corresponds to the correlation enhancement observed in Fig.l. 
The difference of the spin-orbit potential exists mainly in the core region as shown in Fig. 3. 
This is consistent with the charge density distributions shown in Fig. 2. 
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Figure 4: The correlation of the neutron thickness of Pb and the charge radius of ' O. 
ro is the reference radius with ro = 3.0fm. 

Now we investigate the role of the isoscalar-isovector term in the proton-rich oxygen 
isotopes. For the binding energy, the modification due to this term is small. For instance, 
the binding energy of 12 changes by 1.5MeV as the A„ increases from to 0.03. With the 
decrease of the neutron number, the neutron distribution is little changed by the isoscalar- 
isovector term. A correlation between the charge radius of proton-rich oxygen isotopes and 
the neutron thickness of 208 Pb is thus available. The establishment of the data to data 
correlation is helpful since it may bridge the different kind of experiments: the neutron 
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Figure 5: The charge density distribution for 12 O with various cases. 

radius measurement for 208 Pb[10] via parity violating electron scattering and the electron 
scattering on the exotic nuclei[12, 13]. In Fig. 4, the correlation between the charge radius of 
12 ' 13 and the neutron thickness of 208 Pb is shown. The pNN tensor coupling has a much 
less impact on the slope of the correlation relation for 12 ' 13 than that for 23 ' 24 0. As different 
from the neutron-rich isotopes, the correlation relation for these proton-rich isotopes shows 
a nonlinear feature. Contrary to the correlation relation for the neutron-rich isotopes 23 ' 24 0, 
softer models that produce the small neutron thickness for 208 Pb give rise to larger charge 
radii for 12>13 as shown in Fig. 4. 

In Fig.5, the charge density distributions for 12 at different cases are shown. The 
dependence on the isoscalar-isovector term for the charge density distribution is just slightly 
affected by the pNN tensor coupling. Similar to that of the neutron-rich isotopes 23 ' 24 0, the 
change of charge distribution due to the isoscalar-isovector term is concentrated in the core 
region. However, comparing to Fig. 2, the rising or dropping trend of the core charge density 
distribution is different with increasing A„. As shown in Fig. 4, the charge density distribution 
of the drip line isotope ( 12 0) is more sensitive to the change of the density dependence of the 
symmetry energy. This may be attributed to the explicit increasing mean-field expectation 
value of the isovector meson for isotopes towards the proton drip line. The charge density 
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distribution can be related to the cross section of the electron scattering on these exotic 
nuclei. As pointed out in Ref.[14, 22], the form factor in the eikonal approximation at 
high momentum transfer (about 3.0fm _1 ) is sensitive to the charge distribution at the core 
region. According to the calculated charge density distributions in some oxygen isotops, the 
identification of the isoscalar-isovector terms requires the electron scattering on exotic nuclei 
with high momentum transfers. 

We have not made some extensive discussions on the model dependence of the results 
based on various RMF parameter sets instead of using the NL3 parameter set, since a detailed 
discussion on the model dependence may be found elsewhere(for instances, see [3, 7]). For 
the RMF model the variation of the neutron radius of 208 Pb is accomplished by softening the 
symmetry energy through introducing the isoscalar-isovector term that does not change a 
variety of ground-state properties. The different theoretical models may extend the region of 
the uncertainty for the neutron radius of 208 Pb [3, 7] that is about 0.3fm. From the correlation 
between the charge radius of some exotic oxygen isotopes and the neutron thickness of 208 Pb, 
we find there are also considerable uncertainties in the charge radius of exotic oxygen isotopes 
that awaits experimental identification. For instance, even according to an estimation based 
on a linear extrapolation, the uncertainty of the 12 O charge radius that is almost from the 
core region is about 0.3fm. 

In summary, we have studied the charge density distributions of some exotic oxygen 
isotopes in RMF. The isoscalar-isovector term which is not constrained by the present data 
is considered to simulate the uncertainty of the neutron thickness of 208 Pb. The pNN tensor 
coupling that does not affect the symmetry energy in the mean field approximation is included. 
Strong correlations between the charge radius of 12 . 13 . 23 . 24 and the neutron thickness of 208 Pb 
are established, showing that charge radii of 12 . 13 . 23 . 24 o, especially 12 0, exist a considerable 
uncertainty. The neutron thickness for 23 ' 24 reduces with softening the symmetry energy. 
The strong correlation between the radii of 23 ' 24 and the neutron radius of 208 Pb can be 
largely attributed to their special structure and the tensor coupling: The pNN tensor coupling 
explicitly enhances the correlation due to the neutron occupation of the orbital 2si/ 2 - The 
modifications of the charge density distributions for these exotic isotopes due to the isoscalar- 
isovector term are mainly explicit at the core region. The charge radius is sensitive to the 
change of the density dependence of the symmetry energy as the isotope goes towards the 
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proton drip line. These new correlations may bridge two kind of experiments, the neutron 
radius measurement for 208 Pb and the charge distribution measurement for exotic nuclei, to 
jointly constrain the density dependence of the symmetry energy. 
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